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Abstract. Recent UV observations of the globular cluster NGC 2808 (Brown et al. I200H show a significant 
population of hot stars fainter than the zero-age horizontal branch ("blue hook" stars), which cannot be explained 
by canonical stellar evolution. Their results suggest that stars which experience unusually large mass loss on the 
red giant branch and which subsequently undergo the helium core flash while descending the white dwarf cooling 
curve could populate this region. Theory predicts that these "late hot flashers" should show higher temperatures 
than the hottest canonical horizontal branch stars and should have helium- and carbon-rich atmospheres. As a test 
of this late hot flasher scenario, we have obtained and analysed medium resolution spectra of a sample of blue hook 
stars in NGC 2808 to derive their atmospheric parameters. Using the same procedures, we have also re-analyzed 
our earlier spectra of the blue hook stars in ui Cen (Moehler et al. r2UU2l for comparison with the present results 
for NGC 2808. The blue hook stars in these two clusters are both hotter (T c ff >35,000 K) and more helium-rich 
than canonical extreme horizontal branch stars in agreement with the late hot flasher scenario. Moreover, we find 
indications for carbon enhancement in the three most helium-enriched stars in NGC 2808. However, the blue hook 
stars still show some hydrogen in their atmospheres, perhaps indicating that some residual hydrogen survives a 
late hot flash and then later diffuses to the surface during the horizontal branch phase. We note that the presence 
of blue hook stars apparently depends mostly on the total mass of the globular cluster and not so much on its 
horizontal branch morphology. 

Key words. Stars: horizontal branch - Stars: evolution - globular clusters: individual: NGC 2808 - globular 
clusters: individual: NGC 5139 



1. Introduction 



HB (EHB) stars (T off > 20,000 K) having extremely thin 
<0.01M Q ) envelopes. The increase in the bolometric cor- 
rection with increasing temperature turns the blue HB 
into a vertical blue tail in optical colour-magnitude dia- 
grams with the faintest blue tail stars being the hottest 
and least massive (see Fig. ^ for a prominent example of 
a blue tail). It is still unclear how some stars manage to 
lose nearly all of their envelope mass and still undergo the 
helium core flash, as indicated by the long blue tails in 
a number of globular clusters (e.g., NGC 6752, Moehler 
et al.|HP. The globular cluster NGC 2808 represents a 
good template cluster, since it has an extremely long blue 



Low-mass stars burning helium in a core of about 0.5 M 
and hydrogen in a shell populate a roughly horizontal re- 
gion in the colour-magnitude diagrams of globular clus- 
ters, which has earned them the name "horizontal branch" 
(HB) stars. The Galactic globular clusters show a great 
variety in horizontal branch morphology, i.e., in the tem- 
perature distribution of their HB stars. The temperature 
of an HB star depends - at a given metallicity - on the 
mass of its hydrogen envelope, with the hottest or extreme 
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tail extending to very high temperatures and correspond- 
ingly low envelope masses. 

NGC 2808 has been the subject of four recent deep 
ground-based and Hubble Space Telescope (HST) photo- 
metric studies (Sosin et al. 119971 Walker lT999l Bedin et al. 
2000, Brown et al. 1200 1|) . The ground-based studies have 
shown that the blue tail of NGC 2808 extends to extremely 
faint magnitudes (My r* 5™5), with gaps at My « 3 m and 
4™5. These gaps can be seen at V « 18™5 and 20™0 in 
Fig. ^ The brighter gap is thus at a similar position as 
the underpopulated region in the colour-magnitude dia- 
gram of NGC 6752, which separates the EHB stars from 
the classical blue HB stars (Moehler et al. I1997bl 1217001 
Momany et al. 2002;. while the fainter gap coincides with 
the hot end of the blue tail in NGC 6752. The blue tail 
in NGC 2808 thus extends to hotter temperatures than 
that of NGC 6752. That in itself would not be a problem, 
but the spectroscopic analyses of Moehler et al. l|1997bl 
|2"UUU|> show that the blue tail stars in NGC 6752 already 
populate the EHB to the hot end predicted by canonical 
HB models. 




0.5 1 1.5 3 
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Fig. 1. Colour-magnitude diagram of NGC 2808 (Walker 
1999) with our spectroscopic targets marked. Our targets 
were chosen to span the gap at V = 20™0 which, according 
to the hot flasher scenario, separates the canonical EHB 
stars from the fainter blue hook stars. 

Observations of NGC 2808 in the far- and near-UV 
(Brown et al. 120011 their Fig. 3) show that the stars below 
the faint gap form a hook-like feature that extends up to 
0™7 below the hot end of the zero-age HB (ZAHB). Such 
"blue hook" stars are so far known in only three other 
globular clusters: w Cen (D'Cruz et al. 120001 Moehler et 
al. l2UU2|) . NGC 6388 and possibly NGC 6441 (Busso et al. 



2003}, all of which are very massive. Within the frame- 
work of canonical HB theory there is no way to populate 
this region of the UV colour-magnitude diagram without 
requiring an implausibly large decrease in the helium-core 
mass. Thus canonical theory fails to explain both the faint 
UV luminosities and expected high temperatures of the 
blue hook stars. 

Brown et al. |20_01| ) have proposed a scenario to ex- 
plain the blue hook stars. According to this scenario the 
stars brighter than the gap at My ss 4™5 are canonical 
EHB stars which lost all but ^0.01 M Q of their envelope 
mass prior to undergoing the helium core flash. As discov- 
ered by Castellani & Castellani l|1993f) . stars which lose 
more mass than this will leave the red giant branch and 
evolve to high effective temperatures before undergoing 
the helium core flash, producing the so-called "hot flash- 
ers" . Sweigart l|1997|l subsequently showed that hot helium 
flashes can occur either as a star evolves from the tip of 
the red giant branch to the top of the white dwarf cooling 
curve ("early" hot flasher), or later as a star descends the 
white dwarf cooling curve ("late" hot flasher). For even 
higher mass loss, a star will die as a helium white dwarf 
without ever igniting helium. These different evolutionary 
paths are illustrated in Moehler et al. ( 2002 ) and Sweigart 
et al. ^M^ . 

D'Cruz et al. I|1996ll2()()()]l proposed that the blue hook 
stars could be the progeny of such hot flashers, but un- 
fortunately the D'Cruz et al. models were, at most, only 
«0™1 fainter than the canonical ZAHB, much less than 
required by the observations. More recently, Brown et al. 
(2001) have explored the evolution of both the early and 
late hot flashers through the helium core flash to the EHB 
in more detail. Their models show that a late hot helium 
flash on the white dwarf cooling curve will induce sub- 
stantial mixing between the hydrogen envelope and helium 
core, leading to helium-rich EHB stars that are much hot- 
ter than canonical ones, as found previously by Sweigart 
( 1997). This result has been confirmed by the calculations 
of Cassisi et al. (2003), who were able to follow the evo- 
lution of a late hot flasher completely through the helium 
core flash. Brown et al. I|2001[) suggest that such flash mix- 
ing may be the key for understanding the evolutionary 
status of the blue hook stars. Such mixing may also be re- 
sponsible for producing the helium-rich, high gravity field 
sdO stars (Lemke et al. 1997), whose origin is otherwise 
obscure. The flash mixing scenario predicts a gap of about 
6000 K between the canonical EHB stars (i.e., stars with- 
out flash-mixing, including early hot flashers) and the late 
hot flashers as well as a helium dominated atmospheric 
composition for the late hot flashers. The hydrogen defi- 
ciency results in a flux redistribution where more flux is 
emitted shortward of the Lyman edge at 912 A and less 
flux is emitted at ultraviolet (1500 A) and longer wave- 
lengths, which makes flash-mixed stars fainter than canon- 
ical EHB stars at these wavelengths. 

Spectroscopic observations of blue hook stars in lo Cen 
by Moehler et al. ( 2002J showed that these stars reach ef- 
fective temperatures of more than 35,000 K, i.e., well be- 
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yond the hot end of the canonical EHB. In addition, most 
of them show at least solar helium abundances with the 
helium abundance increasing with effective temperature, 
in contrast to canonical EHB stars such as those stud- 
ied in NGC 6752 by Moehler et al. pTO| . Contrary to 
the predictions of Brown et al. ( 2001]) , however, only one 
star showed a helium abundance of log > 0. These re- 
suits may indicate that flash mixing is less efficient than 
assumed, or that some residual hydrogen survives flash 
mixing and then diffuses outward to the surface. It is un- 
clear whether the lower than expected helium abundances 
can still reproduce the UV properties of the stars. UV 
data have been presented for the blue hook stars in the 
central region of ui Cen by D'Cruz et al. i|2(K)0[l . but due 
to crowding these stars are not accessible to ground-based 
spectroscopy. UV data are available for the blue hook stars 
in the outer region of w Cen from the Ultraviolet Imaging 
Telescope (UIT, Stecher et al. I1997fl . but unfortunately, 
these data are too noisy to permit a quantitative test. 

Considering these somewhat ambiguous results and 
the fact that lo Cen is not a typical globular cluster, we de- 
cided to test the late hot flasher scenario by observing the 
more typical globular cluster NGC 2808. We discuss our 
observational data in Sect. and then derive the param- 
eters of the blue hook stars (temperatures, gravities and 
helium abundances) in Sect. [21 In Sect. 0] we compare our 
results with the predictions of the flash-mixing scenario. 

2. Observations and Data Reduction 

2.1. Target Selection 

Our spectroscopic targets were selected from the catalog 
of Walker l|1999l Fig^l, and include seven stars brighter 
than the gap at V « 20 m , and twelve stars fainter than 
the gap. The coordinates and photometry for our targets 
are given in Tabled along with cross- identifications from 
the catalog of Bedin et al. ( 2000 ). 

2.2. UIT observations 

We would like to compare the results of our ground- 
based spectroscopy with the positions of our target stars 
in the UV colour-magnitude diagram. Unfortunately, we 
cannot use the UV colour-magnitude diagram of Brown 
et al. 1)200 1[1 because all of their stars lie within 35" of 
the cluster center and are thus too crowded for ground- 
based spectroscopy. As an alternative, we have derived 
new 1520 A photometry of stars in the outer region of 
NGC 2808 using archival UIT images. The UIT photom- 
etry of NGC 2808 has not been previously published, al- 
though optical spectroscopy of three UV-bright stars on 
the image was presented by Moehler et al. (JT998). The 
UIT image of NGC 2808 was a 979 second exposure ob- 
tained on 19 May 1995 using the Bl filter, which has a 
central wavelength of 1520 A and a width of 350 A. The 
FWHM of the star images was about 5" and the faintest 
hot HB stars were near the sensitivity limit. Therefore, to 



create an ultraviolet colour-magnitude diagram, we per- 
formed 3-pixel (3'.'4) radius circular aperture photometry 
only on isolated stars (all more than 2' from the cluster 
center) which had a positional match with hot (B — V < 
0™2) stars in the catalogs of Walker (^99) or Bedin et al. 
(2OO0J). A star was determined to be isolated if it had no 
neighbours within 5" on either the UIT image, or on the 
lists of hot stars in the optical catalogs. We used the stan- 
dard UIT absolute calibration and determined an aper- 
ture correction of 0™94 using a few isolated UV-bright 
stars. The resulting UIT colour-magnitude diagram ob- 
tained from the 1520 A photometry given in Tabled is 
shown in Fig. [2J and includes all our spectroscopic targets 
except for W26607 which has a blurred PSF on the UIT 
image. The complete ultraviolet photometry and optical 
identifications used in Fig.|3are available upon request to 
the authors. 

Although the UIT colour-magnitude diagram is of 
poorer quality than that shown in Brown et al. (2D0j), it 
is still possible to identify the candidate sub-ZAHB stars. 
Unfortunately, all but one of our blue hook targets are 
less than 0.3 mag fainter than the hottest canonical ZAHB 
model at mi 520 = 16.72. Only W20334 at mi 52 o = 17.35 is 
near the faint end of the blue hook distribution found by 
Brown et al. (2001). Thus our choice of targets is biased 
toward the brighter blue hook stars. 

2.3. Spectroscopy 

We obtained medium-resolution spectra (R «900) of 19 
stars along the blue tail in NGC 2808 with 19 m < V < 
21 m at the VLT-UT1 (Antu) with FORS1 see Fig. EJ). 
The data were obtained in service mode (see Table |5| for 
details). We used the multi-object spectroscopy (MOS) 
mode of FORS1 (slit length 19") with the standard colli- 
mator (O'.'2/pixel), a slit width of 1" and grism B600. As 
FORS1 is equipped with an atmospheric dispersion correc- 
tor, MOS observations at higher airmass are not a prob- 
lem. The detector was a TK2048EB4-1 backside thinned 
CCD with 2048x2048 pixels of (24 ^m) 2 , which was read 
out with high gain (1.46 e _ /count, 5.15 e _ read-out noise) 
and normal read-out speed using one read-out port only. 
This configuration yields a dispersion of 1.2 A/pixel. As 
can be seen from Table |2 the seeing was usually better 
than the slit width, resulting in a seeing-dependent reso- 
lution of the spectra. 

For each night dome flat fields with two different il- 
lumination patterns and CdHeHg wavelength calibration 
spectra were observed. As part of the standard calibration 
we were also provided with masterbias frames for our data. 
The masterbias showed no evidence for hot pixels and was 
smoothed with a 30x30 box filter to keep any possible 
large scale variations, but erase noise. The flat fields were 
averaged for each night and bias-corrected by subtracting 
the smoothed masterbias of that night. From the flat fields 
we determined the limits of the slitlets in spatial direction. 
Each slitlet was extracted and from there on treated like 
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Table 1. Positions and photometric information of target stars 



^2000 


$2000 


Walker (1999) 


Bedin et al. (2000) 


HI152O 








V 


D — V 




1/ 


D — V 










setup 


a 










09 h 12 m 07-20 


-64°54 / 41'.'l 


W6816 


19*35 


+0™04 


B13283 


19*38 


-0*00 


16*51 


09 h 12 m 20?68 


-64°53'09'.'4 


W9301 


19*32 


+0™05 


B12148 


19 m 27 


+0^08 


16*63 


09 h ll m 57 ? 06 


-64°56'18'/6 


W14040 


19*96 


+0™09 


B12955 


19 m 95 


+0^03 


16*75 


09 h ll m 58-81 


-64°56'45'.'6 


W14711 


19*31 


+0™10 


B12583 


19 m 36 


-0*02 


16*48 


09 h 12 m 09-22 


-64° 54' 08'.' 3 


W18803 


20*72 


-0*09 


B20678 


20 m 60 


+0™14 


16*96 


09 h 12 m 09-47 


-64°55'03'.'5 


W18899 


21*00 


-0*07 


B34817 


21 m 06 


-0*08 


17*01 


09 h 12 m 20-95 


-64°53'57'.'2 


W23198 


20*50 


-0*03 


B21073 


20 m 48 


-0*08 


16*71 


09 h 12 m 28-51 


-64°52'19'.'3 


W25794 


20*41 


-0*11 


B21376 


20 m 32 


+0^01 


16*71 


09 h 12 m 31 ? 25 


-64°52'06'.'9 


W26607 


20*58 


+0™03 


B31854 


20 m 71 


+0^00 




setup b 


09 h 12 m 13 ? 42 


-64° 46' 35'.' 4 


W20334 


21™01 


-0*06 


B38604 


21*02 


-0*05 


17*35 


09 h 12 m 17-43 


-64° 48' 03'.' 5 


W21882 


19*75 


-0*05 


B 14334 


19*70 


-0*04 


15*38 


09 h 12 m 29-81 


- 64° 47' 56'.' 1 


W26182 


19*32 


+0™05 


B11516 


19*33 


+0*^03 


16*56 


09 h 12 m 34?ll 


-64°51'23'.'2 


W27412 


19*54 


+0™03 


B15455 


19*64 


-0*00 


16*32 








setup 


c 










09 h ll m 33 ! 72 


-64°56'24'.'9 


W6022 


19*36 


-0*01 


B11345 


19*50 


-0*12 


15*93 


09 h ll m 36-59 


-64°58'01'.'0 


W6849 


20*35 


-0*08 


B23335 


20*35 


-0*13 


16*98 


09 h ll m 37 ? 53 


-64°53'31'.T 


W7084 


20*44 


-0*03 


B26161 


20*59 


-0*03 


16*82 


09 h ll m 42-47 


-64°56'57'.'l 


W8750 


20*62 


-0*10 


B25082 


20*60 


-0*14 


16*87 


09 h ll m 36?34 


-64°54'07'.'7 


W9855 


20*36 


-0*03 


B24578 


20*50 


-0*15 


16*66 


09 h ll m 45-79 


-64°56'27'.'l 


W9863 


20*39 


-0 m 10 


B16326 


20*36 


-0*14 


16*20 



Table 2. Observational parameters 



setup 


start of observation 


seeing 


airmass 


moon 
illumination 


distance 


a 


11/01/2002 05:38:52.722 


0'.'91 


1.319 


0.053 


83-3 




11/01/2002 06:37:32.265 


0'.'61 


1.319 


0.050 


83 ? 4 




15/01/2002 05:13:58.529 


0'.'73 


1.322 


0.027 


95 ? 3 


1) 


20/01/2002 04:50:50.581 


0'.'77 


1.322 


0.350 


105 ? 




08/02/2002 02:26:30.563 


0'.'82 


1.390 


0.161 


85-0 




10/02/2002 02:46:38.787 


0'.'68 


1.354 


0.046 


91?2 


c 


09/02/2002 04:02:45.640 


Of! 67 


1.314 


0.090 


88?2 




09/02/2002 04:46:30.947 


0'.'50 


1.318 


0.088 


88 ? 3 



a long-slit spectrum. The flat fields were normalized with 
4 th or 5 th order polynomials. The dispersion relation was 
obtained from the wavelength calibration frames by fitting 
3 rd to 4 th order polynomials to the line positions along the 
dispersion axis. We used 13 to 16 unblended lines between 
3600 A and 6200 A and achieved an r.m.s. error of 0.05 A 
to 0.08 A per row. 

Due to exposure times of 2700 sec the scientific obser- 
vations contained a large number of cosmic ray hits. Those 
were corrected with the algorithm described in Gossl & 
Riffeser piOljl . using a threshold of 15cr and a FWHM of 
1.5 pixels for the cosmic ray hits. As the routine is not orig- 
inally intended for the use with spectra, we also reduced 
the uncorrected frames to allow a check for any possi- 
ble artifacts of the cosmic ray cleaning procedure. The 
slitlets with the stellar spectra were extracted in the same 
way as the flat field and wavelength calibration slitlets. 
The smoothed masterbias was subtracted, and the spectra 



were divided by the corresponding normalized flat fields, 
before they were rebinned 2-dimensionally to constant 
wavelength steps. We then filtered the uncleaned frames 
along the spatial axis with a median filter of 7 pixels width 
to erase cosmic ray hits. We identified regions uncontami- 
nated by any stellar source and approximated the spatial 
distribution of the sky background by a constant. For the 
uncleaned frames the sky background was fit on the me- 
dian filtered frames and subtracted from the unfiltered 
ones, whereas for the cleaned frames both sky fit and sub- 
traction were performed on the unfiltered rebinned data. 
The sky-subtracted spectra were extracted using Home's 
(1986) algorithm as implemented in MIDAS. Finally the 
spectra were corrected for atmospheric extinction using 
the extinction coefficients for La Silla (Tug ll97?j) as im- 
plemented in MIDAS, because they provide the closest 
approximation to Paranal conditions, for which no spec- 
troscopic extinction coefficients are available. 
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Table 3. Atmospheric parameters of the stars in 
NGC 2808 



Fig. 2. Ultraviolet- visual colour-magnitude diagram of 
NGC 2808 as derived from UIT observations. The spec- 
troscopic targets are marked by big dots. The zero- age 
HB (solid line) is computed for [M/H] = —1.1 and has 
been transformed assuming a distance modulus (m-M)o 
of 15™04 and a reddening of Eb-v = 0™18 as suggested 
by Bedin et al. l|2T)0Ti|l for the Zinn & West metallicity 
scale. If we adopt the distance and reddening determined 
by Bedin et al. (2000) for the Caretta & Gratton scale, 
then the ZAHB would be ss0™2 fainter. 



For a relative flux calibration we used response curves 
derived from spectra of LTT 3864, LTT 4364, and 
LTT 6248 with the data of Hamuy et al. l(TM?|) . The 
response curves were fit by splines and averaged for all 
nights. The individual target spectra were corrected for 
any Doppler shifts determined from Balmer and/or helium 
absorption lines. These shifts varied more than expected 
from one observation to the next. This can be understood 
by the fact that the seeing disk was smaller than the slit 
width, so that small shifts in the position of the stars 
within the slit can introduce shifts in the wavelength. A 
shift of O'.'l corresponds to a shift of 0.6 A. Any remain- 
ing cosmic ray hits in the stellar spectra were corrected 
by comparing the individual spectra of the stars. Due to 
the seeing being smaller than the slit width spectra of the 
same star from different nights can have different resolu- 
tion. In order to get a well defined resolution of the co- 
added spectra, we convolved spectra obtained with better 
seeing to the resolution determined by the worst resolu- 
tion of each setup. This resulted in resolutions of 5.4 A 
(0"9) for setup a, 4.8 A (Of! 8) for setup b and 4.0 A (Of! 7) 
for setup c. 



Star 


■J 

X 




T e ff 

[K] 


logg 


log 

o n H 


W6022 


2.23 


30300± 


480 


5.57±0.05 


-1.99±0.07 


W6816 


2.39 


21300± 


540 


4.88±0.07 


-2.17±0.05 


W6849 


1.88 


37500± 


330 


5.77±0.09 


-0.82±0.05 


W7084 


1.63 


35300± 


580 


5.91±0.10 


-0.57±0.05 


W8750 


2.52 


48000±1300 


6.21±0.14 


-1.78±0.15 


W9301 


2.14 


21100± 


560 


4.92±0.07 


-2.05±0.05 


W9855 


2.45 


35700± 


580 


5.60±0.12 


-0.74±0.05 


W9863 


2.04 


53000±2100 


5.81±0.10 


-1.56±0.14 


W14040 


2.16 


40400± 


220 


5.63±0.05 


-0.31±0.03 


W14711 


2.92 


20100± 


520 


4.87±0.07 


-1.93±0.05 


W18803 


1.70 


35700± 


510 


5.96±0.09 


-0.38±0.03 


W18899 


1.91 


35500± 


790 


5.70±0.12 


-0.21±0.05 


W20334 


2.15 


36300± 


900 


5.79±0.15 


+1.02±0.02 


W21882 


4.09 


67000±1500 


6.99±0.09 


-2.91±0.12 


W23198 


2.38 


34600± 


400 


5.81±0.07 


-1.22±0.05 


W25794 


1.44 


36600± 


460 


5.67±0.09 


-0.65±0.03 


W26182 


1.79 


24200± 


640 


5.22±0.07 


~2.37±0.05 


W26607 


2.91 


36300± 


520 


5.63±0.09 


-0.36±0.03 


W27412 


2.64 


28700± 


640 


5.37±0.07 


-2.77±0.12 



Table 4. Atmospheric parameters of the stars in u> Cen 



Star 


x 2 


T ff 


logg 


t> n H 






[K] 




WF3-1 


1.34 


36000±2100 


5.97±0.40 


-2.07±0.50 


BC6022 


1.96 


46000±1100 


6.14±0.15 


-1.76±0.14 


BC8117 


1.65 


30000±1600 


5.43±0.21 


-2.33±0.26 


BC21840 


2.94 


36000±1000 


5.44±0.19 


-0.78±0.12 


C521 


2.20 


35000± 620 


5.86±0.12 


-0.88±0.07 


D4985 


2.50 


38000± 700 


5.79±0.17 


-0.80±0.14 


D10123 


2.40 


35200± 690 


5.73±0.14 


-0.84±0.09 


D10763 


2.85 


36000±1700 


4.82±0.15 


+0.50±0.03 


D12564 


1.75 


36800± 960 


5.73±0.17 


-0.50±0.05 


D 14695 


2.16 


38800± 230 


5.31±0.12 


-0.16±0.05 


D15116 


1.50 


40000±1700 


6.08±0.38 


-0.51±0.17 


D16003 


1.45 


36600± 550 


5.89±0.12 


-1.00±0.09 



3. Analysis 

The spectra, given in Fig. [21 show a large variety of helium 
line strengths. Part of this variation is due to variations in 
effective temperature, but - as we will see in the analysis 
- part is also due to variations in the atmospheric helium 
abundance along the blue tail in NGC 2808. 

3.1. Atmospheric parameters 

We fitted the hot stars (T cff >33,000 K) with H-Hc non- 
LTE model atmospheres to derive effective temperatures, 
surface gravities, and helium abundances. The helium- 
rich non-LTE model atmospheres were calculated with a 
modified version of the accelerated lambda iteration code 

mm 



of Werner & Dreizler (1999). The model atoms for hy- 
drogen and helium as well as the handling of the line 
broadening for the spectrum synthesis are similar to those 
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Fig. 3. Spectra of blue tail stars in NGC 2808 sorted by their number in the photometry of Walker (1999 ). The spectra 
were normalized by eye to allow a better comparison. The lines mark the position of the He I and He n lines. 



of Werner (ljT99*6t) . The calculation of the helium-poor version 1 of the LINFOR program (developed originally by 

non-LTE model atmospheres is described in Napiwotzki Holweger, Steffen, and Steenbock at Kiel University). To 

(1997). For the cooler stars we used ATLAS9 model atmo- establish the best fit, we used the routines developed by 

spheres for solar metallicity (Kurucz ll993fl to account for Bergeron et al. I|1992f) and Saffer et al. ( 1993), as modified 

effects of radiative levitation (see Moehler et al. [2000 for by Napiwotzki et al. (1999), which employ a x 2 test. The 

details), from which we calculated spectra with Lemke's — _ 

For a description see http://a400.sternwarte.uni-erlangen.de/~ai26/linfit/ 
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a necessary for the calculation of \ 2 is estimated from 
the noise in the continuum regions of the spectra. The 
fit program normalizes model spectra and observed spec- 
tra using the same points for the continuum definition. 
During the analysis of the data presented here we real- 
ized that helium-poor and helium-rich spectra required 
different sets of continuum points. Especially the use of 
the continuum points derived from helium-poor spectra 
for the analysis of helium-rich spectra can introduce large 
systematic errors, esp. overestimates of the helium abun- 
dance due to continuum points in the wings of strong he- 
lium lines and/or too narrow fitting windows for helium 
lines. We therefore refined the definition of the continuum 
points to correct for such errors. Using this new defini- 
tion of the continuum points, we re-analysed the data for 
the blue hook stars in uj Centauri presented in Moehler et 
al. (2002). The new results differ from the old ones usu- 
ally within the mutual error bars, except for the two most 
helium-rich stars, as expected. 

We used the Balmer lines H/3 to H§ (excluding H £ to 
avoid the Can H line), the Hei lines AA 4026 A, 4388 A, 
4471 A, 4921 A and the Hen lines AA 4542 A, 4686 A for 
the helium-poor stars. For the helium- rich stars we also 
included the Hei lines AA 4713 A, 5015 A and 5044 A in 
the fit. The results are given in Table El for NGC 2808 and 
Table H for uj Cen. 




7.6 r i i i i i i i i i i i i i i i i i i i ~ 

75000 70000 65000 60000 
T eff M 

Fig. 4. Comparison of the stellar parameters for W21882 
(solid square) with evolutionary tracks for helium (dashed 
curves) and carbon-oxygen (solid curves) white dwarfs. 
The masses of the tracks are 0.423, 0.448 and 0.480 M© for 
the helium white dwarfs and 0.491, 0.494, and 0.529 M Q 
for the carbon-oxygen white dwarfs in order of increasing 
gravity. 

As can be seen from Table [3] W2 1882 shows both a 
much higher effective temperature (67,000 K) and surface 
gravity (logg = 6.99) than all other stars. These stellar pa- 
rameters indicate that W21882 may be a low-mass white 
dwarf, just in the mass range between the most massive He 



white dwarfs and the least massive C/O white dwarfs (cf. 
FigEJ). This suggests an absolute magnitude of My «6™7, 
which would place W21882 at a distance of at most 4 
kpc (zero reddening assumed), compared to 9.6 kpc for 
NGC 2808. Unfortunately, the radial velocity estimates 
derived from our spectra are too crude to determine clus- 
ter membership. However, as can be seen from Table|21 the 
model fit to W21882 has by far the largest \ 2 value of all 
our targets. We examined the three individual spectra to 
search for variability that might indicate binarity or prob- 
lems with the spectra, but did not find any such evidence. 
Additional observations of this very hot star are needed to 
better determine its cluster membership and evolutionary 
status and we therefore omit it from all further discussion. 
Figs. and El show the results for all other stars. 
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Fig. 5. Atmospheric parameters derived from the spectra 
of stars along the blue tail/blue hook in NGC 2808 (filled 
squares, Table |3| compared to HB evolutionary tracks. 
Also shown are blue tail stars from NGC 6752 (open tri- 
angles, Moehler et al. 2000| and the blue hook stars in uj 
Cen (open squares, Table 0J. The tracks for an early hot 
flasher (long-dashed line) and a late hot flasher (short- 
dashed line) show the evolution of such stars from the 
zero-age HB (ZAHB) towards helium exhaustion in the 
core (terminal-age HB = TAHB). The solid lines mark 
the canonical HB locus for [M/H] = —1.5 from Sweigart 
( 1997). The dotted line connects the series of ZAHB mod- 
els computed by adding a hydrogen-rich layer to the sur- 
face of the ZAHB model of the late hot flasher. The large 
dots mark - with decreasing temperature - hydrogen layer 
masses of 0, 10~ 7 , 10~ 6 , 10~ 5 , 10~ 4 M Q . 

In order to verify if the atmospheric parameters can 
reproduce the UIT observations, we calculated UIT fluxes 
from the observed V magnitudes, using the model spectra 
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Fig. 6. Helium abundances vs. effective temperature for 
the stars along the blue tail/blue hook in NGC 2808 (filled 
squares, Table 0J. Also shown are blue tail stars from 
NGC 6752 (open triangles, Moehler et al. |2000) and the 
blue hook stars in uj Cen (open squares, Table H}. 



from the best fits. The model spectra were first reddened 
by Eb-v = 0™18. Then mean fluxes through both the 
UIT 1520 A and the V filters were calculated to obtain 
the mi52o — V colour, which was added to the observed 
V magnitude from Walker (1999) to obtain the expected 
TO1520 magnitude. The results are compared to the UIT 
observations in Fig. [7| The agreement is good for the ma- 
jority of the stars considering the low quality of the UIT 
data. The largest discrepancies are for the two helium-rich 
stars (W6849 and W14040), which have much fainter ob- 
served UIT magnitudes than predicted. A possible cause 
for these discrepancies might be the omission of carbon 
or other metal opacity in our model atmospheres, which 
Lanz et al. |2003j have shown can affect both the derived 
atmospheric parameters from optical spectra, and the ul- 
traviolet flux distribution. 

3.2. Carbon abundances 

Along with a large helium enhancement, the flash-mixing 
scenario also predicts a strong surface carbon enhance- 
ment. The flash-mixed model computed by Cassisi et al. 
(2003) is strongly enhanced in both carbon (mass frac- 
tion 0.029) and nitrogen (mass fraction 0.007). The corre- 
sponding solar mass fractions are 0.0028 and 0.00081 for 
C and N, respectively. We see lines of C in AA 4070, 4170, 
4650 A in several of our spectra, but our low resolution 
and S/N prevents a detailed abundance analysis. 

To obtain a rough estimate of the carbon abundance, 
we began with the H-He NLTE models used to derive Teff, 




±±i 1520, cale. 

Fig. 7. Observed UIT magnitudes vs. UIT magnitudes 
calculated from observed V magnitudes and model spec- 
tra. The solid line marks identity between observed and 
calculated UIT magnitudes. Filled squares mark helium- 
rich stars, open squares mark helium-poor stars. 

log g and the helium abundance. Keeping these parame- 
ters fixed we then computed five additional models with 
CNO 2 logarithmic mass fractions relative to solar of —1.5, 
-0.5, 0.0, +0.5, and +1.0. As noted by Lanz et al. pTOljl . 
this approach is inconsistent for large carbon abundances, 
in that the additional C opacity can alter the atmospheric 
structure significantly, and therefore the derived stellar 
parameters. The impact of this additional C opacity has 
so far been explored only in the case of the field He-sdB 
star PG1544+488. For this star Lanz et al. 1|2003[1 found 
that the value of logg given by NLTE models with a C 
abundance of 0.02 by mass, as obtained by fitting the UV 
C III lines, was about 0.5 dex larger than the value given 
by models without C. Additional calculations are needed 
to determine the carbon abundance at which it is possible 
to neglect carbon in the model atmosphere and still de- 
rive accurate stellar parameters. We note, however, that 
we find a supersolar carbon abundance in only one star 
(W20334), and so our crude approach is likely adequate 
for the remaining stars. 

The comparison between model spectra and observa- 
tions resulted in a carbon abundance compatible with the 
cluster abundance ([Fe/H] = —1.15, Harris IT!)96[) except 
for the three most helium-rich stars W14040, W18899, 
and W20334. For W14040 (log^= -0.31) and W18899 
(log 2 ^ 1 = —0.21) we estimate logarithmic carbon mass 
fraction relative to the Sun of —0.5 and 0.0, respectively. 
The most helium-rich star, W20334 (log +1.02), 

2 Although oxygen should not be enhanced in flash-mixed 
stars, we do not expect the enhanced oxygen in our models to 
alter our estimate of the carbon abundance. 
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shows Cm lines at AA 4070 A and 4650 A, which are 
even stronger than predicted by a model atmosphere with 
a CNO abundance of ten times solar. In this case, the 
limitations of our two-step approach (see above) become 
evident. As noted above, a more accurate abundance de- 
termination for this star requires an analysis with a model 
grid taking carbon opacity into account for the determi- 
nation of T e ff and logg (Lanz et al.[2003). 

A similar estimate of the carbon abundances for the u> 
Cen stars was not possible due to the poor S/N and reso- 
lution of these data. Nevertheless, the strong C in features 
found in some of the u> Cen spectra suggest an enhanced 
carbon abundance in some of the blue hook stars. 



4. Discussion 

Our analysis of blue tail and blue hook stars in NGC 2808 
shows that these stars form two different groups accord- 
ing to their helium abundance (cf. Fig. |8}: All stars 
brighter than the gap at V«20 m are helium-poor, although 
they cover a large temperature range (cf. Fig. |HJ. The 
stars fainter than the gap at V«20 m , on the other hand, 
show mostly super-solar abundances, with three excep- 
tions: W8750, W9863, and W23198. Fig. ©shows that the 
helium-rich stars cluster between 35,000 K and 40,000 K, 
well beyond the hot end of the canonical EHB. We will 
now discuss the evolutionary status of the helium-poor 
and helium-rich stars in more detail. 



1 - 



.-1 



20 



21 



Fig. 8. Helium abundance vs. V magnitudes. The gap in 
the colour-magnitude diagram is located at V « 20 m , and 
log £^ = — 1 corresponds to solar helium abundance. 



4.1. Helium-Poor Stars 

As can be seen from Fig. [31 the six helium-poor stars 
cooler than 31,000 K occupy a range in effective temper- 
ature and surface gravity similar to the blue tail stars in 
NGC 6752. Their helium abundances are also similar to 
those observed in NGC 6752, indicating that the gravita- 
tional settling of helium has been equally efficient in both 
clusters. Thus, the helium-poor stars brighter than the 
gap at V«20 m agree with the predictions of canonical HB 
evolution. 

The two helium-poor stars hotter than 40,000 K, 
W8750 and W9863, appear to be in a post-HB evolution- 
ary state This can be seen in Fig. |5| where we compare 
the stellar parameters of these stars with HB and post- 
HB evolutionary tracks for canonical models, flash-mixed 
models with a thin hydrogen surface layer, and flash-mixed 
models with no envelope hydrogen. W9863 appears to lie 
along the hottest canonical post-HB track at a point where 
the evolution is relatively slow. However, things are less 
simple for W8750. While it lies close to the post-HB track 
of a late hot flasher with a 10 -5 M Q hydrogen surface 
layer, it also lies at a point where the evolution is rapid. 
Moreover, it is difficult to understand why W8750 has a 
helium-poor surface when the other stars with possible 
flash mixing are much more helium-rich. In order to place 
it on the least massive canonical track (where evolution is 
slow in this temperature range), its gravity would have to 
be overestimated by about 0.4 dex, which is about 3o\ 

W23198 is a less clear case: In Figs. Eland it lies at 
the cool end of the population of helium-enriched stars, 
but well separated from the canonical EHB stars. Judging 
from its position in Fig. 9, W23198 may be a late hot 
flasher with very high hydrogen content and/or effective 
diffusion or an example for shallow mixing (see below). 

4.2. Helium-Rich Stars 

The parameters of the helium-rich stars provide support 
to the late hot flasher scenario: Figs. [5] and show that 
a gap between «31,000 K to «35,000 K separates the 
helium-rich from the helium-poor stars, as predicted by 
the flash-mixing scenario. The HB track for the early hot 
flasher in Fig. passes through the temperature gap, but 
the evolution is then very fast, making a contamination of 
the temperature gap by such stars unlikely (Moehler et al. 
2002). In addition, Fig.|f)|shows the expected clustering of 
the helium-rich stars between 35,000 K and 40,000 K. 

Contrary to the predictions of Brown et al. (2001) 
and Cassisi et al. ('2003, but consistent with our previ- 
ous results for oj Cen, the atmospheres of the blue hook 
stars still show some hydrogen. This result has been dis- 
cussed by Cassisi et al. I|200 3 ) , who find that to reproduce 
even the highest observed helium abundances, they have 
to reduce the efficiency of the flash mixing by a factor 
of about 20,000. However, this reduction only applies if 
the observed helium abundances reflect the actual helium 
abundances in the envelopes of the blue hook stars. 
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Fig. 9. Comparison of the stellar parameters for W8750, 
W9863, and W23198 (solid squares) with EHB and post- 
EHB evolutionary tracks. Each track is represented by a 
series of points separated by a time interval of 5 TO 5 yr in 
order to indicate the rate of evolution. The solid line in the 
upper righthand corner denotes the canonical ZAHB. The 
two lowest gravity tracks are for canonical, i.e., unmixed 
models with hydrogen-rich envelopes. The two interme- 
diate tracks labeled H) are for flash-mixed models with 
a hydrogen-rich surface layer containing either 10~ 5 M Q 
(higher gravity track) or 10~ 4 M Q (lower gravity track). 
The highest gravity track labeled He-C) is for a flash- 
mixed model with no envelope hydrogen. 



This apparent discrepancy could be explained if some 
residual hydrogen survived flash mixing and later diffused 
outward, thus producing a thin hydrogen-rich layer at the 
surfaces of the blue hook stars. Such diffusive processes 
arc believed to be responsible for the low helium abun- 
dances of the sdB stars and are estimated to operate on 
a time scale much shorter than the HB lifetime. It is un- 
clear, however, if diffusion will turn a helium-rich star into 
a helium-poor one. Groth et al. (19831 found that atmo- 
spheric convection (which would work against diffusion) 
can exist in hot subdwarfs if the helium abundance is suf- 
ficiently high. Also mass loss might affect the atmospheric 
abundances of the blue hook stars by reducing the effi- 
ciency of the hydrogen diffusion. The low helium abun- 
dances observed for canonical EHB stars, however, sug- 
gest that mass loss does not prevent the formation of a 
hydrogen surface layer. The range in the hydrogen abun- 
dances of the blue hook stars might indicate that varying 
amounts of hydrogen survive flash mixing or that the ef- 
ficiency of diffusion differs from star to star. In any case 
the high helium abundances observed in some of the blue 
hook stars would be difficult to understand if their at- 
mospheres were not enriched in helium during the helium 
core flash. The increase in the mean atmospheric helium 
abundance with increasing effective temperature is also 
consistent with flash mixing. As can be seen in Fig. El 



adding a hydrogen layer to the surface of a late hot flasher 
moves the track towards lower temperatures and gravities 
(see Moehler et al. 2002 for more details). The addition of 
a hydrogen layer of < 10~ 4 M Q would actually improve 
the agreement between the predicted and observed tem- 
peratures of the blue hook stars while at the same time 
preserving the temperature gap between these stars and 
the canonical EHB stars. 

Another possibility for explaining the hydrogen abun- 
dances of the blue hook stars has been discussed by Lanz 
et al. ( 2003), who found that there are two types of flash 
mixing: "deep" and "shallow" , depending on how far the 
envelope hydrogen is mixed into the core during a late 
helium-core flash. During deep mixing the envelope hy- 
drogen is mixed all the way into the high temperature 
regions near the site of the helium flash and therefore 
rapidly burned. This is the type of mixing discussed by 
Brown et al. (200TJ and Cassisi et al. |2003). A blue hook 
star following deep mixing will be helium- and carbon- 
rich with very little hydrogen left in its atmosphere. In 
contrast, during shallow mixing the envelope hydrogen is 
only mixed with the outer layers of the core. Since the 
temperatures within these layers are too low for proton- 
capture nucleosynthesis, all of the envelope hydrogen will 
survive the flash-mixing phase. The surface composition of 
a blue hook star following shallow mixing will therefore re- 
main hydrogen-rich, although its surface composition will 
be diluted by helium- and carbon-rich material from the 
core. Unfortunately the calculations of Lanz et al. (2003) 
indicate that shallow mixing only occurs over a very nar- 
row range of mass loss in metal-poor stars. Thus, while 
shallow mixing may be a viable explanation for the resid- 
ual hydrogen found in some field He-rich sdB stars (Lanz 
et al. I2003f) . it is probably of minor importance for the 
blue hook stars in NGC 2808. A detailed study of metal 
abundances might be able to distinguish the signatures of 
shallow mixing from those of diffusion. 

5. Conclusions 

The high temperatures and high helium abundances re- 
ported here for the blue hook stars in NGC 2808 and 
to Cen provide general support for the flash-mixing sce- 
nario of Brown et al. (12001[) . The three most helium-rich 
stars in NGC 2808 also show evidence for carbon enrich- 
ment. However, all of our targets show some hydrogen in 
their atmospheres, and only one target in each cluster has 
log 0. This is partially a selection effect, as our tar- 
gets are preferentially the brighter blue hook stars, and the 
most helium-rich star in each cluster is either the faintest 
(W20334 in NGC 2808) or second-faintest (D10763 in lu 
Cen). However, it is also likely that some hydrogen sur- 
vives a late hot flash, and subsequently diffuses to the 
surface during the HB phase. 

Moehler et al. (1997a) reported an isolated helium- 
rich sdB star in M15 with the expected properties (T e ff 
= 36,000 K, M v = +4 1 ?7, log^= +0.82) of a flash- 
mixed star. However, significant numbers of blue hook 
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stars have been found only in the most massive Galactic 
globular clusters: u Cen (My = -10™29), NGC 2808 
(My = -9™39), NGC 6388 (My = -9'°42), and pos- 
sibly NGC 6441 (My = -9*64), all data for My from 
Harris (1996). Other globular cluster in the same mass 
regime are 47 Tuc (My = -9*42), NGC 2419 (My = 
-9*58), and NGC 6715 (My = -10 m 01). The optical 
colour-magnitude diagrams of NGC 2419 (Harris et al. 
1997) and NGC 6715 (Momany, priv. comm., Rosenberg 
et al. l2003[) show stars faint enough to be blue hook stars. 
The situation for 47 Tuc is unclear due to contamination 
by the Small Magellanic Cloud. 

Apparently the number of blue hook stars depends on 
the total mass of the globular cluster and not on the num- 
ber of other hot HB stars (see Rosenberg et al. 20Q3|for 
an independent discussion of this effect). For instance, 
NGC 6752, which is a factor of 100 less luminous than 
NGC 2808 and has a very long and well populated blue 
tail, contains 59 EHB stars with 3 m < My < 4™5 (pho- 
tometry of Buonanno et al. I1986|) . compared to 34 such 
stars observed in NGC 2808 (photometry from Walker 
1999). Even more EHB stars can be seen in the data 
of Momany et al. QMty for NGC 6752. However, while 
NGC 2808 has 28 blue hook stars with 4 m 5 < My < 5™5, 
no such stars are found in NGC 6752. Likewise, the UV 
colour-magnitude diagram of Brown et al. |20D1| shows 
that the population of blue hook stars dominates over the 
population of EHB stars in NGC 2808 (46 vs. 29 stars, 
respectively). As the mass loss ranges populating the hot 
end of the canonical extreme HB and the blue hook region 
are not too different, it is hard to understand why globu- 
lar clusters with large numbers of canonical EHB stars do 
not show at least some blue hook stars as well. 
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